A Dynamical Mean-Field Study of Ferromagnetism 
in the Periodic Anderson Model 
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The ferromagnetic phase diagram of the periodic Anderson model is calculated using dynamical 
mean-field theory in combination with the modified perturbation theory. Concentrating on the in- 
termediate valence regime, the phase boundaries are established as function of the total electron 
density, the position of the atomic level and the hybridization strength. The main contribution to 
the magnetic moment stems from the /-electrons. The conduction band polarization is, depending 
on the system parameters either parallel or antiparallel to the /-magnetization. By investigating the 
densities of states, one observes that the change of sign of the conduction band polarization is closely 
connected to the hybridization gap, which is only apparent in the case of almost complete polariza- 
tion of the /-electrons. Finite-temperature calculations are also performed, the Curie temperature 
as function of electron density and /-level position are determined. In the intermediate-valence 
regime, the phase transitions are found to be of second order. 
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I. INTRODUCTION 

The periodic Anderson model (PAM) represents proba- 
bly the simplest way of modeling some gfnthe rich physics 
found in Lanthanides and ActinidestlB'El. One believes 
that most of the physical properties typical for these ma- 
terials originate from the interplay between the incom- 
pletely filled 4/ or 5/ shells which contain almost local- 
ized electrons, and a broad conduction band of s, p and 
d-electrons. In the periodic Anderson model, this is simu- 
lated by an atomic-like level with strong on-site Coulomb 
interaction which hybridizes with an uncorrelated con- 
duction band. Usually the model is further simplified by 
assuming both the atomic- and conduction states to be 
non-degenerate. The on-site Coulomb interaction within 
the /-states has to be considered as the largest energy 
scale of the system since the electrons occupying these 
states are less screened than the conduction electrons. In 
the PAM, the Coulomb interaction is further necessary 
to ensure that the /-levels are not fully occupied, thus re- 
sembling the situation found in the Lanthanides and Ac- 
tinides. Due to the incompletely filled /-levels, local mo- 
ments could be formed at every lattice site. It has been 
a long-standing question whether these moments order 
magnetically or whether the local moments are.-screened 
by conduction band electrons (Kondo screening)!! In this 
paper we will focus on those parameters where the / elec- 
tron density is non- integer (intermediate valence regime) . 
In this parameter regime, we find a stable ferromagnetic 
phase and investigate its properties. 

Kondo screenir^Jj^a n |^cii^ | h£ | _subject of exten- 
sive investigatipnJraiBmffl'ril. Besides few ex- 
act statementst3li3, several approximate or numerical 
methods have been applied to explore the possibil- 
ity of ferromagnetism in the PAM, e. g. Hartree-jEock 
calculationsO, spectrajL,£te»sity approach (SDA)E3o, 
slave-boson teckniauesE3Eirn3 and dynamic mean-field 
theory (DMFT)Lrcil The first methods have severe lim- 
itations concerning the description of the low-energy 
physics. The DMFT-based method promise an improve- 



ment at this point. However, there are still only few 
results available on ferromagnetism in the PAM within 
dynamical mean-field theory. 

The antiferromagnetic- .phase -of— tiic PAM has been 
discussed in more detaiflLjEjO'ESSo. It seems to be 
widely accepted that close to the symmetric parameter 
conditions, i. e. half-filling and symmetric DOS, antifer- 
romagnetism has to be expected. This can also be con- 
cluded from the following: ThepPAM can be mapped 
onto an effective Hubbard modelliSiij. In the vicinity of 
the symmetric point the effective Hubbard model will 
also be close to half-filling. In that case the Hubbard 
model is-commonly expected to show antiferromagnetic 
orderingcJ, therefore the PAM will also have a tendency 
towards antiferromagnetic ordering as was indeed foundE. 
Ferromagnetism is expected further away from this sym- 
metric point. A ferromagnetic phase was found for the 
PAM both in the KondojcawQtaBB and in the in- 
termediate valence regionlI§II3c3. Whereas in the former 
case, the driving force towards ferromagnetism can be 
understood via an effective Heisenberg-like coupling of 
the /-spins (RKKY interaction)!], such a picture is not 
easily applied to the IV regime due to the non-integer 
density of /-electrons and their effective itineracy. 

In this paper, we will investigate ferrorti agneti sm in 
the PAM using dynamical mean-field theoryt^E3Ej. This 
theory is based on the assumption of a /c-independent 
self-energy, which becomes exact in the limit of infinite 
dimensions (d — oo)c3. As pointed out in reference |2^, 
in this limit the lattice coherence and the exchange ef- 
fects due to the Pauli principle are preserved, contrary 
to other approaches basec -QfL.jj expansion (N being the 
degeneracy of the model)E3td. Especially the exchange 
effects should be vital for ferromagnetism, therefore the 
DMFT seems a method of choice. 

In the next section, the DMFT together with the mod- 
ified perturbation theory (MPT) are introduced. The 
DMFT will lead to a mapping of the PAM onto a single- 
impurity model, which then will be solved by applying 
the MPT. The results concerning zero temperature as 
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well as finite temperatures will be presented and dis- 
cussed in section III. 



summary here. Starting point is the following ansatz for 
the self-energyaCJ: 



II. 



DMFT AND THE MODIFIED 
PERTURBATION THEORY 



V a {E) = U(n [ H) + 



(SOC) 



(4) 



The periodic Anderson model is defined by its Hamil- 
tonian 

h =£ <k) S y^ + E e fflfi° + (!) 

i.a i.a 

(fie-) and st (fj a ) are the annihilation and creation 
operators for an electron in a non-degenerate conduc- 
tion band state (localized /-state), and n^y — fj a fia is 
the occupation number operator for the /-states. The 
dispersion e(k) describes the propagation of free, i. e. 
unhybridizcd conduction electrons, e/ is the position of 
the free /-level relative to the center of mass of the con- 
duction band density of states. The hybridization V is 
taken as a real, /c-independent constant, and finally U is 
the Coulomb repulsion between two /-electrons on the 
same lattice site. 

The quantity of interest will be the /-electron Green 
function 
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a a and /3 CT are introduced as parameters to be deter- 
mined later. T,^° C \e) is the second-order contribu- 
tion to jiciJjiii'bation theory around the Hartree-Fock 
solutior£j'El'c3. The Hartree-Fock solution introduces 
another free parameter, namely the chemical potential 
within this calculation: ft. A priori it is not clear whether 
this should be equal to the chemical potential in the full 
(DMFT-MPT) calculation, or whether e. g. it should be 
determined such that the electron density on the impu- 
rity site of the SIAM is equal both for the Hartree-Fock 
and DMET-MPT calculation. In reference |3^ and other 
papersr 2 lr 3 l yet another condition aias used to determine 
ft. There, the Lwttinger theorema or equivalently the 
Friedel sum ruleOcj was forced to hold. Since these the- 
orems are applicable only for T = 0, this limits the cal- 
culations to zero temperature. In order to access finite 
temperatures, we used the condition of identical electron 
densities for the Hartree-Fock and the full calculation 
(ni/' H = Tie). With this choice, the Friedel sum rule 
is still fulfilled in a large parameter region as could be 
shown for the pure SIAMEj. A more detailed analysis 
of the different possibilities to determine ft is found in 
reference 

Next, the remaining parameters a a and f3 a have to be 
determined. Instead of using the "atomic" limit of V = 
•as done e. g. in references [L2,p6l§^ we make use of 



(2) 



joi-cmploy the dynamical 
It was shown that 



To determine this functior 
mean-field theory (DMFT)^I 
in the limit d — > oo, the self-ene rgy Eg (E) becomes 

purely local, i. e. A-independentoEj. In this case, the 
self-energy of the PAM is equivalent to the self-energy 
of a prpcedy. defined single-impurity Anderson model 
(SIAM)E3'BEl. The latter has to be defined by the so- 
called self-consistency condition 



^tnc'momcnts of the spectral density 

M<»> = [dEE n A(f\E) = ([[...[f a ,HU...,HUft} + ) 



-fold commutator 



(5) 



A^\E) = --^G ( il{E + l Q + ) 



A a (E) = E - ( e/ - /i) - £ ff (£) - ( G\H (E) 



instead of the usual definition A = V ? 

pure SIAM. Using 
in references 34 3 



(3) 



for the 



lerturbation theory it has been shown 
that for a three-dimensional system, 
the local approximation (equivalent to the limit d = oo) 
provides already for satisfactory results. 

Now one is left with the problem to solve the SIAM 
defined by equation (||). Here we use the modified per- 
turbation theory (MPT) which can be understood as 
an improvement of the IPT scheme introduced in ref- 
erences |2^|^_Ihis method was presented in more de- 
tail elsewhereHll'tS, so we will restrict ourselves to a short 



where [^4;B]_ ([A;£?] + ) denotes the commutator (anti- 
commutator) of the operators A and B and 3a; denotes 
the imaginary part of x. As indicated, these can be cal- 
culated on two different ways, therefore conditions (sum 
rules) can be derived. To determine a a and (3 ai the first 
four moments n € {0, ... 3} have to be used since the 
n = 0- and n = 1-moments are fulfilled for any a a and 
f3 a . In the n = 3- moment a higher-order correlation func- 
tion that we will call bandshift B a is introduced: 



( n (/))(l-( n (/)))(B CT -e / ) = 

= E v ^t/^ (/ ]-i)> 



(G) 



= -ha J dE f-(E)A<j(E) - 1) g\2(e) 

with the Fermi function f-(E) — (exp(fSE) + 
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Now, the solution of the SIAM by the MPTiuis to be 
integrated into the DMFT-self-consistency loopz3: Start- 
ing with a guessed value of T, a (E) , equations (]|) and (||) 
are evaluated, and then the new MPT self-energy (|J) is 
calculated for the appropriate SIAM. This procedure is 
iterated until a self-consistent solution is found. Within 
this formalism, the spontaneous symmetry-breaking nec- 
essary for ferromagnetic solutions can be introduced via 
spin-asymmetric starting values for the self-consistency 
cycle. In general, it is not necessary to introduce a 
symmetry-breaking zero-field. 

At this point, let us comment on the quality of our ap- 
proximative method: Being based on perturbation the- 
ory, the MPT is expected to give reliable results for small 
interaction strengths. For the symmetric SIAM, it was 
even showncj that the perturbative expansion is essep, 
tially equivalent to the exact Bethe-ansatz solutionc§c3 
and therefore also valid in the strong coupling regime. 
For asymmetric parameters, we improve on the pertur- 
bation theory by enforcing the correct high-energy ex- 
pansion for the self-energy and equivalently the Green 
functions. This is archieved by fulfilling the first four 
sum rules defined by the spectal moments (||) which leads 
to the correct determination of the spin-dependent posi- 
tions and weights of the charge excisions up the or- 
der jj in accordance with reference |50pJ. From this, we 
derive our proposition that the MPT can also be rea- 
sonably applied to a PAM in the intermediate-to-strong 
coupling regime. Although for the high-energy (high- 
temperature) behaviour this follows from the discussion 
above, the quality of the special low-energy properties of 
the PAM within the MPT is not known a priori. To es- 
timate the significance of the results for the intermediate 
coupling strengths, comparison with exact or numerically 
exact methods is necessary. For example, one can apply 
the MPT to a SIAM without the context of the DMFT 
and compare with the exactly known properties of that 
modelu. This comparison was done in reference ^ and 
the results can be summarized as follows: The charge ex- 
citations are at the proper positions for symmetric as well 
as asymmetic parameters. Concerning the low-energy be- 
haviour, there is a qualitative agreement but quantitative 
discrepancies: The Friedel sum rule is fulfilled in a large 
parameter space off the symmetric point, but the Kondo 
temperature does not follow an exponential law. A power 
law is found instead. To summarize, the MPT does qual- 
itatively include essential parts of the Kondo physics, but 
is prone to deviations concerning energy (temperature) 
scales. Further tests of the MPT can be done by com- 
paring results for the PAM with different methods based 
on DMFT, aSj-fen g. the numerically exact quantum Monte 
Carlo (QMC)Btl, exact diagonalization (ED)E3 or numeri- 
cal renormalization group theory (NRG)Ej. Comparisons 
of this kind for theparamagnetic PAM have been pub- 
lished in reference [jj] for QMC and reference |52] for ED 
calculations. Also the results of reference [l4| can easily be 
compared with above-mentioned references. Again, the 
qualitative features compare well, however temperature 



scales concerning low-energy properties cannot be repro- 
duced quantitatively. Also, the ^-dependence seems rto 
be overestimated by the MPT when compared to NRGE3. 
We do believe that these shortcomings of the MPT do not 
inhibit the analysis of the ferromagnetic properties of the 
PAM. In the discussion of our results, we will substanti- 
ate this claim. 

In this paper, we use the following system parame- 
ters: The conduction band is described by a semi-elliptic 
density of states of unit width centered at E = 0, thus 
defining the energy scale. The position of the /-level, 
e / will be given relatively to the center of gravity of the 
conduction band. The hybridization is taken to be con- 
stant. Unless otherwise noted, we chose V = 0.2. The 
Coulomb interaction strength U is set to U = 4 for most 
calculations, thus representing the largest energy in the 
system. The total number of electrons and the 

temperature T will be varied. The latter will be given 
in K/eV within the energy scale defined by the width of 
the conduction band. 



III. RESULTS AND DISCUSSION 

Before discussing the results concerning ferromag- 
netism in the PAM, let us recall some results for the 
paramagnetic case as they will be important for the fol- 
lowing. The density of states (DOS) consists of the two 
charge excitation peaks approximately at e/ and e/ + U. 
In addition, a sharp resonance close to the chemical po- 
tential is induced by the correlations, the Kondo reso- 
nance. Due to lattice coherence, the Kondo resonance is 
split by the coherence gap in case of symmetric param- 
eters. Another result for the paramagnetic PAM found 
within the DMFT-MPT scheme is the existence-of local 
Kondo singlets for large hybridization strengths^. More 
detailed reviews on the paramagnetic PAM can be found 
in references 0g||Jpgg|| 

In figure ^, the region of ferromagnetic order is plotted 
in the e/ vs. n^ *-* phase diagram. We have restricted our 
evaluation to the so-called intermediate valence region 
(IV region) by positioning the /-level within the conduc- 
tion band which extends from —0.5 to 0.5. This results in 
a /-occupancy n^> smaller than unity. Previous investi- 
gations have show n . thatj ferromagnetism can exist in this 
parameter regimet§c3E3. The opposite case is the so- 
called Kondo regime, which is obtained by taking e/ well 
below the lower conduction band edge, and U sufficiently 
large so that ef + U is clearly above the conduction band. 
This leads to a nearly half-filled /-level. As discussed in 
the introduction, a half-filled /-level should introduce a 
tendency towards antiferromagnetism which we have not 
yet investigated. In the following we will therefore con- 
centrate on fcrromagnctism in the intermediate valence 
region, but we will also be noting some trends on how 
the ferromagnetic phase would continue into the Kondo 
regime. As can be seen in figure [|, the region of ferro- 
magnetic order is quite large, determined by an upper 
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bound for e/ and upper and lower bounds for n^ tot \ As 
will be shown below, the ferromagnetic region continues 
into the Kondo regime. 

In the left panel of figure 0, the /- (s-) magnetization, 

n (/ ' s) -n (/ ' s) 

m^' s > — If „ s — jtjnrr is plotted as function of /-level po- 
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sition e/. The parameters are as in figure |lj, the band 
occupation is given by n^ tot ^ = 1.3. Ferromagnetic or- 
der breaks down for e f approaching the center of grav- 
ity of the conduction band as it is commonly expectedE3. 
Within our numerical accuracy, this quantum phase tran- 
sition is of second order. 

In the other limit, i. e. e/ below the lower edge of the 
conduction band (e/ < —0.5), the ferromagnetic solution 
remains stable. We confirmed that down to ef = —1.2 
the magnetization behaves as one would extrapolate from 
figure p|. So, the ferromagnetic phase plotted in figure |l| 
continues down into the Kondo regime and connects to 
the ferromagnetic phase found in reference [?]. There, the 
authors found a stable ferromagnetic phase for n( tot ) < 
1.6. 

The right picture of figure || shows the magnetization 
as function of the hybridization strength for e/ = —0.4, 
the remaining parameters as in the left panel. The ferro- 
magnetic order is destroyed by large hybridizations. This 
can be understood in terms of local singlet formation. 
If the conduction electrons screen the magnetic moment 



"i — i — r 



"i — i — i — r 



"i — i — i — r 




j L 



of the localized electrons, as discussed in detail in refer- 
ence [l4], magnetic ordering will npt .be possible any more. 
Comparison with NRG results£3E3 seem to suggest that 
the ^-dependence of certain "low-energy" quantities is 
over-estimated by the MPT, as e. g. the size of the co- 
herence gap. This implies the possibility that the critical 
V is underestimated since the formation of local Kondo 
singlets is due to the low-energy physics of the PAM. 
For V — > the ferromagnetic phase is stable, for nu- 
merical reasons no calculations for V < 0.03 have been 
performed. For V — 0, the /-level and the conduction 
band decouple completely. The /-level basically corre- 
sponds to the zero-bandwidth Hubbard modelEU in which 
ferromagnetism is not stable. In figure Q, the thin lines 
represent the conduction band polarization . In gen- 
eral, it is found to be antiparallel to However, as 
will be shown below, a parallel alignment of m^- 1 and 
m,( s ) is also possible for certain parameters. 

In figure ^, the (/-dependence is examined. The inset 
shows that ferromagnetism sets in at a critical U c and 
saturates quickly afterwards. The absolute value of U c 
varies strongly with ef. On shifting e/ below the con- 
duction band into the Kondo regime, U c seems to sat- 
urate at a value of Uc K « 2.5. However, there is no 
true saturation, because when e/ is low enough so that 
the upper charge excitation is close to or even below the 
chemical potential, the physics of the system will change 
again. On increasing e/, U c increases as well and exhibits 
a sharp upturn at e/ — > —0.2. In reference f^, a ferromag- 
netic phase was found in the Kondo regime of the PAM 
for interaction strengths smaller than our . As for 

the hybridization, the MPT seems to misjudge the en- 
ergy scale. We will argue below that the Kondo scale is 
decisive for ferromagnetism in the Kondo regimd!, but 
not in the intermediate valence regime. 




FIG. 1: e f vs. n (tot) phase diagram for T = 0, U = 4, 
V = 0.2 



FIG. 2: /- (s) -magnetization as function of the hybridization 
strength V and /-level position e/ for T = 0, n^ ^ = 1.3 and 
(7 = 4 plotted in thick (thin) lines (left panel: V — 0.2, right 
panel: e/ — —0.4). 
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-0.5 -0.4 -0.3 -0.2 

e f 

FIG. 3: Critical interaction strength U c as function of ef 
for T = 0, V = 0.2 and n (tot) = 1.3. The inset shows 
the respective /-magnetization as function of U for e/ = 
-0.5, -0.3, -0.25 (from top to bottom). 
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FIG. 4: /- and s-magnetization as function of electron density 
for various e/ and U = 4, V = 0.2 and T — 0. 



region depend one/. The lower boundary is in fact de- 
termined by the number of /-electrons which itself 
depends on e/ and n,( tot ). The phase transition occurs 
for each et at that specific electron density, where 
drops below 0.56 ± 0.02. Approaching the half-filled sys- 
tem, no ferromagnetism is found. This is similar to the 
findings in reference 0. However as disc^sa^Lahove, an- 
tiferromagnetic order is to be expectecloLlcj'Ej'cfl in this 
case. 

For ef < —0.32 the ferromagnetic phase can be sep- 
arated into two regions of different properties. In the 
low-density regime, the /-polarization tends to satura- 
tion, the s-polarization is positive. On the contrary, 
in the high-density region (1.1 < n( tot ) < 1.7), the /- 
polarization is almost constant, below saturation, and 
the conduction band polarization is of opposite sign (an- 
tiferromagnetic coupling between /- and s-band). For 
more clarity, we will call the first case parallel, the latter 
case antiparallel coupling of /- and s-electrons. This lat- 
ter case of antiparallel coupling is generally expected for 
the PAM due to its close relation to the Koado, lattice 
model via the Schrieffer- Wolff transformation^^. This 
transformation maps the PAM onto a model with an an- 
tiferromagnetic interband-spin exchange. One condition 
for the Schrieffer- Wolff transformation is that the /-level 
should be half-filled, which is not met here. However, 
for large total occupation number, rv-f' is closer to unity 
than for small n( tot ). Therefore the finding of parallel 
coupling for low, and antiparallel coupling for high elec- 
tron densities is no contradiction to the Schrieffer- Wolff 
transformation. 

In figures || and |[ the densities of states are plotted 
for various band occupations with e/ = —0.4, V = 0.2, 



U = 4 and zero temperature. The respective left panel 
corresponds to a projection onto /-, the right panel onto 
s-states. The upper charge excitation, located approxi- 
mately at E « et + U = 3.6 is not visible in these figures. 
The arrows indicate the respective position of the chem- 
ical potential. To get more insight into the complicated 
structure of the DOS, we start by explaining figure |^, 
where we have forced the system to be paramagnetic. In 
figure ^, the stable ferromagnetic DOS are shown and 
will be explained below. 

Starting with the /-DOS in figure |[ one observes three 
structures: the first is the lower charge excitation located 
roughly at the lower band edge. With increasing elec- 
tron density, this peak shifts to lower energies. The next 
structure is the Kondo resonance which is in the vicin- 
ity of the chemical potential ji. Due to the asymmetry 
of the parameters, i. e. the two charge excitation peaks 
are not symmetric around the chemical potential and the 
band occupation is well below n( tot ) = 2, the Kondo res- 
onance is not exactly at E = fi but slightLy above. This 
is in agreement with previouSr£ndingsE3o and was also 
observed in theJ>IRG methocO and the single-impurity 
Anderson modeltlH. With increased electron density, the 
Kondo resonance moves towards higher energies, this is 
clearly related to the according shift of the chemical po- 
tential [i. It is reasonable to expect the Kondo reso- 
nance to vanish for decreasing n^' . Within the MPT, the 
Kondo resonance remains stable for all parameters dis- 
cussed in this paper. Only for much reduced where 
et is located well above /i, the Kondo resonance disap- 
pears. As already mentioned, there is the "coherence 
gap" in the center of the Kondo resonance for symmetric 
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FIG. 5: Paramagnetic /- and s-densities of states for ej — 
-0.4, V = 0.2, U = 4 and T = 0. The different panels 
correspond to different electron densities n' tot ' 
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FIG. 6: same as figure ^, but showing the (stable} ferromag- 
netic solution. The parameters are as in figure H, thus cor- 
responding to the dotted line of figure ^| Solid line: spin |; 
dotted line: spin J.. 



parameters. Due to the asymmetry of the parameters 
in figure |, this gap is closed. Only for the n^ ^ = 1.4 
case (upper panel of figure ||) , there is a slight dip in the 
/-DOS visible. This complies with the findings of ref- 
erence [l3| The third feature in the /-DOS is a broad 
structure representing states induced by the hybridiza- 
tion with the conduction band. A fourth feature is of 
course also present in the DOS, although not plotted in 
figures | and | the upper charge excitation which is lo- 
cated roughly at e f + U = 3.6 and has a FWHM (full 
width at half maximum) of 0.05. 

All these structures have their counterpart in the 
conduction-band DOS shown in the right panel of fig- 
ure ^| The charge excitation induces via the hybridiza- 
tion states in the s-DOS. Where the Kondo resonance is 
visible in the /-DOS, the s-DOS shows also some anoma- 
lies. The original conduction band can also be clearly 
recognized. 

Let us now discuss how the DOS are modified in the 
ferromagnetic phase, fn figure ^ the (stable) ferromag- 
netic DOS are plotted for the same parameters as in fig- 
ure |5|. The solid lines represent the spin-f, the dotted 
lines the spin- J, DOS. For the lowest panel (n' tot ^ = 0.6), 
the system is paramagnetic. For the ferromagnetic so- 



lutions, the following picture emerges: Firstly, the lower 
charge excitation is fully polarized, i. e. it consists purely 
of spin-| states. Correspondingly, the upper charge exci- 
tation, which is not visible, exists only in the spin- j DOS. 
Some interesting changes occur to the Kondo resonance: 
Whereas in the spin-| channel it disappears almost com- 
pletely, there is still a peaked structure visible in the 
spin- J, channel. Contrary to what one would expect for 
the Kondo resonance, its shift towards higher energies 
on increasing n( tot ) is stronger than the according shift 
of [i. It is therefore not located at /z any more. Although 
this structure clearly develops from the Kondo resonance, 
its physical interpretation (Kondo screening) cannot so 
easily be transfcrcd. Finally there are some substantial 
changes in the s-DOS as well: Due to the full polar- 
ization of the lower charge excitation, the hybridization- 
induced s-states in the same energy range are also purely 
of majority spin. The feature that corresponds to the 
Kondo resonance is, for similar reasons, only visible in 
the spin-j-s-DOS. Finally, in the spin-| DOS, a new gap 
appears which separates the lower charge excitation from 
the remainder of the bare band. Due to this new gap, 
the system becomes a semi-metal around quarter-filling 
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( n (tot) _ i meanm g that the spin-f-DOS vanishes at 
H, so only spin-j-electrons contribute to electrical cur- 
rent. The new gap, which we call hybridization gap is 
not equivalent to the previously discussed coherence gap 
although the latter is of course also induced by the hy- 
bridization. We will explain the differences between these 
gaps below. This gap is also present in the /-DOS, but 
not visible in the plotted figures because of an artificial 
broadening of the lower charge excitation which is neces- 
sary for numerical reasons. 

From these observations the change of sign of the con- 
duction band polarization can be understood: In the case 
of low electron density, the chemical potential lies within 
the lower charge excitation. As explained above, due to 
the hybridization there are only spm-| conduction band 
states available in the occupied region of the DOS, the 
conduction band polarization is necessarily parallel to the 
/-polarization. But on the other hand, with increasing 
electron density, fi is increasing as well. The chemical 
potential now has to cross the hybridization gap which 
exists only in the spin-f DOS. As a result, there are more 
occupied spin-j than spin-| states and the conduction 
band polarization is antiparallel to the /-magnetization. 
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FIG. 7: Temperature dependence of the DOS for n (tot) = 1.2. 
The remaining parameters are as in figure a. 



Let us now discuss the difference between the well- 
known coherence gap and the above-described hybridiza- 
tion gap. The coherence gap is found for a symmetric 



PAM with on-site hybridizations. Contrary to the single- 
impurity case, a gap develops in the center of the Kondo 
resonance located at E = /i. Shifting away from sym- 
metric parameters, this coherence gap gets suppressed, 



Honance is still visible 
So why is the new 



although a dip within the K 
for a large parameter regio 

gap observed in figure |^ not this coherence gap? This 
question is answered by the fact that the new gap does 
not shift with n on increasing the electron density. Al- 
though there is no proper Kondo resonance in the case 
of ferromagnetic order, those structures that clearly de- 
velop from it do shift with increasing n^ tot \ Within this 
structure, even a dip which is the relic of the coherence 
gap is visible for n( tot ) = 1.4. The position of the new 
gap is neither related to /z nor the relics of the Kondo res- 
onance, it is rather located at e/ . And what is the origin 
of the new gap? Since its position is determined by the 
value of e/ and its width scales with V, it reminds o 
gap found e. g. within the SDAEa or alloy analog; 
but also for the C/_p case of the PAM which is ex- 
actly diagonalizabletl. This indicates that this gap is due 
to level-repulsion between the lower /-peak and the con- 
duction band, as the coherence gap can be interpreted 
as level-repulsion between a virtual /-level at the fermi 
energy and the conduction band. To distinguish between 
these gaps, we call the former "hybridization" and the 
latter "coherence" gap. 

At first sight, it surprises that this gap was pat fa 
any paramagnetic DMFT-based calculatiorJlMJ-t 
This is in our opinion simply due to damping effects 
caused by the correlation, i. e. a finite imaginary part 
of the self-energy. In case of saturated ferromagnetism, 
vanishes and the corresponding spin-f DOS 
should resemble the U = DOS including the hybridiza- 
tion gap at e/. In the spin- J.- or also the paramagnetic 
DOS, all structures are broadened due to the finite damp- 
ing. This effectively closes the hybridization gap for these 
DOS. This proposal can be supported by investigating 
the temperature-dependence of the DOS. On increasing 
the temperature, the magnetization will decrease. The 
system is pushed away from saturation, damping effects 
occur also for spin-f electrons. In figure ^, the DOS for 
n (tot) _ 2 2 are plotted for various temperatures (the re- 
maining parameters as in figure |J). For these parameters, 
the Curie temperature is T c w 408. One can see how the 
gap closes on approaching T c . More complicated is the 
case of n < - tot - ) = 0.8 in figure ^. On the one hand, ferro- 
magnetic saturation is not reached (cf. figure [|) and the 
hybridization gap is partly closed by damping. And on 
the other hand, the chemical potential is located in the 
same energy range, further complicating the situation. 

In the following we will discuss the temperature de- 
pendence in more detail. Figure || shows magnetiza- 
tion curves for n( tot ) = 1.2 and three different values 
of e/ € {-0.5,-0.4,-0.3} (from top to bottom). The 
lowest line (e/ = —0.3) shows a clear second-order phase 
transition. The middle case corresponds to the data of 
figure The phase transition at T c w 408 is within 
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FIG. 8: Magnetization curves for U = 4, V = 0.2, 
and various e/ £ {—0.5, —0.4, —0.3} (from top to bottom) 



numerical accuracy of second order. The demagnetiza- 
tion process is dominated by a transfer of spectral weight 
from the minority- to the majority-spin spectrum of the 
lower charge excitation peak as can be seen in figure |7]. 
However, for the upper charge excitation, plotted in fig- 
ure ^| for the same parameters as figure [?], not only the 
reverse transfer of spectral weight occurs but also a spin- 
dependent bandshift can be observed. When decreasing 
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FIG. 9: Upper charge excitation peak of the /-DOS for the 
same parameters as in figure |?] 



e/ towards the lower band edge of the conduction band, 
the phase transitions become sharper, and finally, first 
order phase transitions are found for < —0.5 within 
numerical accuracy. This is similar to findings for the 



Hubbard model0, where phase transitions were found to 
be of second order for lower, and of first order for higher 
electron densities. In the PAM, the density of correlated 
electrons, increases on lowering e/, therefore the 

behaviour is analogous to the findings for the Hubbard 
model. As discussed in reference [si], we cannot exclude 
that this is an artifact of our method. In figure 0, the 



500 




FIG. 10: Curie temperatures as function of e/ for U 
V = 0.2. The electron density is as indicated. 



4 and 



Curie temperature is plotted as function of the /-level 
position. The trends are very similar to the T = 0- 
magnetization as seen in figure 0. For e/ approaching 
the upper critical value, both T c and m(T = 0) vanish. 
In the opposite direction, both quantities keep increasing, 
however, with diminishing slope. The Curie temperature 
as function of band occupation is shown in figure [fl] for 
e/ = —0.4 and e/ = —0.3. 

As already mentioned in the introduction, the mag- 
netism in the Kondo limit of the PAM can be,piiderstood 
in terms of an effective RKKY interactionLrEil between 
the /-spins. However, in the intermediate-valence regime 
under investigation here, the effective RKKY exchange 
should be suppressed due to the non-integer /-electron 
filling and the effective itineracy of the /-electrons as 
indicated by the broader /-peaks in the DOS. We be- 
lieve that instead of the interband-RKKY exchange, a 
"narrow-band" mechanism is the main driving force to- 
wards ferromagnetism in the IV regime. The effective 
itineracy of the /-electrons due to the hybridization in 
combination with the strong on-site Coulomb interac- 
tion leads to t he .| Stro ,ng similarities to the well-known 
Hubbard modelEBE^CJ. Recently the existence oLfetm- 
magnetism was shown within the Hubbard modelOOLj. 
One important condition for it is the shape of the free, 
i. e. U = DOS entering the model. A highly asym- 
metric DOS with a divergence or sharp peak close to 
the baad edge enhances the possibility of ferromagnetic 
orderEH. Turning back to the periodic Anderson model we 
note that this condition favoring ferromagnetic order is 
excellently met for the hybridization-broadened /-level 
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in combination with the hybridization-induced /-states 
within the conduction band. Our proposal, that ferro- 
magnetism in the IV regime of the PAM originates from 
an intraband mechanism is supported by several findings: 

i) Whereas in the Kondo regime, small values of U can 
already lead to a ferromagnetic phase, in the IV regime, 
significantly larger values of U are necessary. 

ii) The particular low-energy physics, that are decisive for 
ferromagnetism in the Kondo regime, have no significant 
influence on the magnetism in the IV regime. 

iii) The critical n^ tot \ where ferromagnetism breaks 
down, is in fact determined by a critical number of cor- 
related /-electrons, which is of similar magnitude as for 
the Hubbard model. 

iv) The polarization of the conduction band seems to 
have almost no influence on the magnetic properties in 
the IV regime. 

Let us discuss these points in more detail: The first 
point follows directly from figure ^ and reference ^[ In 
the Kondo regime, where the. low-energy physics deter- 
mine the magnetic propertiesQ, the MPT probably over- 
estimates the absolute value of U c . However, the sharp 
increase of U c for increasing e/ indicates the impor- 
tance of a strong-coupling mechanism in the IV regime. 
The second point can be made clear by comparing the 
MPT results with those obtained by two other methods, 
namely the spectral density approximation (SDA)li3 and 
the modified alloy analogy (MAA)Ej. Whereas all three 
methods do fulfill the first four moment-sum r ules (|) 
important for the strong-coupling limit, the MPT surely 
represents a major improvement over the other two meth- 
ods concerning the low-energy physics. The SDA and 
the MAA completely neglect the Kondo physics from the 
very beginning while the MPT includes, as discussed in 
section II, at least qualitatively the special low-energy 
properties of the PAM. Now, if the low-energy physics 
have a major influence on ferromagnetism in the IV 
regime of the PAM, one would expect a dramatic change 
of the properties of the ferromagnetic phase. However, 
in all three methods, the ferromagnetic phase turns out 
to be very similar, key features as e. g. the phase di- 
agram, but also the unusual conduction band polariza- 
tion are present in all approximations. This clearly in- 
dicates the minor importance of the special low-energy 
physics of the PAM when describing ferromagnetism in 
the IV regime. The contrary seems to be the case for 
the Kondo regime: here-the SDA does not give a sta- 
ble ferromagnetic phasdl§E3, whereas within the MPT, 
the ferromagnetic solution remains stable in the Kondo 
regime and the parameters leading to ferromagnetism 
agree very well with the results of reference if one ig- 
nores the already discussed overestimation of U c . The 
third point can be clarified by pointing out that for an 
RKKY mechanism, the density of the conduction band 
electrons should have a major influence on the effective 
exchange integrals, and therefore also on phase bound- 
aries. However, as discussed in the context of figure ^, 
the lower phase boundary as function of 

n (tot) Jg 

governed 



by the number of /-electrons and is independent of . 
This is the expected behaviour for an intraband mecha- 
nism as proposed by us. It is clear that for the Kondo 
regime with fixed n^' ~ I, the situation is different. Fi- 
nally let us discuss the influence of the conduction band 
polarization. It is clear that it is decisive for an RKKY 
mechanism. However, our results indicate that the mag- 
nitude and the unusual behaviour of mS s ' have no direct 
effect on the stability of ferromagnetism in the IV regime 
(cf. figures |, | and[n]). This is a clear indication that 
an RKKY mechanism is suppressed here. 

To summarize the discussion, we believe that two dif- 
ferent mechanism are to be considered when discussing 
ferromagnetism in the PAM: In the Kondo regime, an ef- 
fective RKKY exchange (interband mechanism) is mainly 
responsible for ferromagnetism, in the IV regime, the 
RKKY exchange is suppressed and an intraband mecha- 
nism similar to the one found for ferromagnetism in the 
single-band Hubbard model should be of major impor- 
tance. 

From this proposal it follows that the bandshift (|^), 
which is introduced by the fourth moment of the spec- 
tral density, should also be decisive for ferromagnetism 
in the PAM in theJV regime as it was shown to be in 
the Hubbard modeEI. The low-energy physics, however, 
are vital in the Kondo regime since the RKKY exchange 
is believed to origin in the formation of Kondo screening 
cloudso but less important in the IV regime. Therefore 
the shortcomings of the MPT in reproducing the low- 
energy scale of the PAM are nonrelevant for the present 
investigation. 

Of course, both the RKKY and the intraband mech- 
anism are always present, only each respective impor- 
tance varies. Whereas in the Kondo limit, the effective 
/-electron itineracy becomes negligible, in the IV regime, 
strong-coupling "Hubbard-like /-band" picture is domi- 
nant. This explains the smooth transitions found be- 
tween the to scenarios (cf . figures [| || and [To]) . 




FIG. 11: Curie temperatures as function of electron densifr 
for different e / . The remaining parameters are as in figure fLl 
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IV. SUMMARY 

In this paper we have presented a detailed investigation 
of ferromagnetism in the intermediate valence regime of 
the periodic Anderson model. We applied the dynamical 
mean-field theory in combination with the modified per- 
turbation theory. We have shortly discussed the quality 
of the MPT based on a comparison with exact results for 
the SIAM and numerically exact results for the paramag- 
netic PAM. Although the MPT cannot recover the cor- 
rect exponential energy scale for the low-energy physics, 
the qualitative features of the low-energy physics emerge 
properly. And the high-energy features are reproduced 
with much higher credibility. 

We have established the phase diagram as function 
of the total electron density and the position of the /- 
level. Furthermore, an upper bound of the hybridiza- 
tion strength was found, above which no magnetic or- 
dering seems to be possible due to local Kondo singlet 
format ior£3. The critical interaction strength U c was 
found to vary strongly when moving from the intermedi- 
ate valence to the Kondo regime. The main contribution 
to the magnetic moment originates from the localized 
/-electrons. However, the conduction band is, due to 
the hybridization, also polarized. As function of electron 
density, the conduction band polarization changes sign. 
For higher electron densities, the conduction-band polar- 
ization is antiparallel to the /-level magnetization mS^ 
as one would expect from the Schrieffer- Wolff transfor- 



mation. For low densities, this transformation is not ap- 
plicable. Here, the conduction band polarization is par- 
allel to The densities of states in the ferromagnetic 
state show rather complicated structures. Besides relics 
of the Kondo resonance and the dip corresponding to the 
coherence gap, another gap is present which we named 
hybridization gap. This gap is closed when approaching 
the paramagnetic state. Investigating the temperature 
dependence, we find phase transitions of second and first 
order. We also have presented the Curie temperatures as 
function of /-level position and total electron density. 

Concerning ferromagnetism in the Kondo and Inter- 
mediate valence regime, two different pictures emerge: 
Whereas in the Kondo regime, an Kondo-screening in- 
duced RKKY exchange leads to the ordering of the £.- 
spins already for relatively low interactions strengthsQ, 
in the intermediate-valence region, another mechanism 
becomes more important: due to the effective itineracy 
of the /-electrons, the /-levels now represent an strongly 
correlated narrow band, which by itself can lead to fer- 
romagnetism as known for the Hubbard model. 
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